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Abstract. We provide a theoretical model that describes the dielectric coupling of a 2D layer 
of graphene, represented by a polarization function in the Random Phase Approximation, 
' and a semi-infinite 3D substrate, represented by a surface response function in a non-local 

■ formulation. We concentrate on the role of the dynamic response of the substrate for low- 

c/3 | frequency excitations of the combined graphene-substrate system, which give rise to the 

stopping force on slowly moving charges above graphene. A comparison of the dielectric 
loss function with experimental HREELS data for graphene on a SiC substrate is used to 
estimate the damping rate in graphene and to reveal the importance of phonon excitations in 
an insulating substrate. A signature of the hybridization between graphene's ir plasmon and 
the substrate's phonon is found in the stopping force. A friction coefficient that is calculated 
' O ' for slow charges moving above graphene on a metallic substrate shows an interplay between 

^ • the low-energy single-particle excitations in both systems. 
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1. Introduction 

After the initial wave of intense studies of electronic and transport properties of graphene 
[1, 2], the focus of research seems to be turning to graphene's interactions with a dielectric 
environment [2, 3]. While it has been recognized early on that the presence of charged 
impurities in the material surrounding graphene is a likely cause of its peculiar minimum 
conductivity [4], much attention has also been paid to the dielectric screening of these 
impurities by either their host material or various high-kappa dielectrics [5]. Besides an 
obvious interest in Si0 2 as the most common substrate for graphene that is produced by the 
exfoliating technique [6, 7], interactions of epitaxially-grown graphene on a SiC substrate 
have also attracted considerable attention [8, 9]. Most recently, several reports have appeared 
studying the properties of graphene on metallic substrates [10, 11], as well as graphene's 
interaction with organic solvents and ionic solutions [12, 13]. 

In view of these developments and possible extensions of the recent Electron Energy 
Loss Spectroscopy (EELS) studies of free-standing graphene structures [14, 15] that would 
include the presence of a substrate, it is of interest to study the interaction of graphene with 
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externally-moving charged particles in a dielectric environment. Although it was recently 
shown that a dielectric environment can exert strong effects on the plasmon excitations of 
carbon nanotubes due to coupling of their a + n and ix plasmons with high frequency modes 
in the nearby material [16], in this paper we limit ourselves to the effects of a substrate 
on slow-moving external charges such that only the low-energy excitations of graphene's 
7r electron bands occur. From an applied point of view, this problem is highly relevant to 
a recent experiment that used a High-Resolution EELS (HREELS) technique to probe the 
momentum-resolved low-frequency plasmon excitations of doped graphene grown on a SiC 
substrate [17]. Moreover, interactions with slow heavy particles are relevant to studies of 
intercalation of alkali-metal atoms [18], friction forces on migrating atoms and molecules 
moving near graphene [19, 20], and monitoring of ion flow in an aqueous solution adjacent to 
carbon nano-structures [21]. 

From a theoretical point of view, the interaction of graphene with slow external 
charges can be described by the recently developed dielectric function in the Random Phase 
Approximation (RPA), in which graphene's ix electron bands are treated in the approximation 
of linearized electron energy dispersion, giving rise to the picture of massless Dirac fermions 
(MDF) in two dimensions (2D) [22, 23]. This theoretical model captures all of the physical 
processes relevant to such interactions, including graphene's 2D 7r-plasmon mode and both 
the intra-band and inter-band single particle excitations (SPEs). Furthermore, the MDF-RPA 
dielectric function for graphene takes into account the strong dependence of these excitation 
modes on the level of doping in graphene, which often results from charge transfers between 
graphene and the substrate [6, 17, 11] or may be induced by a gate potential. However, all 
RPA dielectric functions suffer two major shortcomings: they neglect the local-field effects 
(LFE) due to electron correlations, and they have an inherent assumption of an infinite lifetime 
of electron excitations. Only one of these shortcomings can be qualitatively corrected in 
the MDF-RPA dielectric function for graphene at a time, e.g., by using either the Hubbard 
approximation (HA) for the LFE [2] or by introducing a finite decay rate, 7, using Mermin's 
procedure [24, 25]. 

With respect to graphene's coupling with a substrate, we note that few theoretical studies 
have gone beyond treating the substrate in a static mode and most assume a vanishing 
graphene-substrate gap, with the exception of recent studies of the remote scattering of charge 
carriers in graphene on surface phonon modes of a strongly polar substrate [26, 27, 28]. 
Therefore, in this work we first provide a detailed description of graphene's coupling with 
a substrate by combining the LFE or Mermin-corrected MDF-RPA dielectric function for 
graphene with the surface response function for a substrate, which is described by a bulk 
dielectric function that is allowed to depend on both the frequency and the three-dimensional 
(3D) wavenumber. We note that the dependence of the substrate's dielectric function on the 
wavenumber is often neglected, giving rise to the local approximation for its surface response 
function that is appropriate for describing cases such as the non-dispersing Fuchs-Kliever 
surface phonons in strongly polar substrates [27, 28]. However, when the substrate's bulk 
dielectric function depends strongly on the 3D wavenumber, such as in the case of a metal, 
one has to resort to the so-called Specular Reflection Model (SRM) [29] in order to obtain 
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a surface response function that takes into account the spatial dispersion of the substrate's 
excitation modes, and hence gives rise to a non-local description of the substrate. 

While such a formulation of the graphene- substrate interaction allows for the possibility 
of studying the hybridization of collective modes in the two systems [16], we focus on 
analyzing the effects of the graphene- substrate dielectric coupling on the energy losses of a 
slow external charge moving parallel to graphene. It has been experimentally confirmed that 
slow ions can lose their kinetic energy during grazing scattering from an insulating surface 
due to excitations of optical phonons in the substrate [30], and it would therefore be interesting 
to explore how such a process would change in the presence of graphene and the possibility 
of simultaneous excitations of graphene's n plasmon and its intra-band SPEs. Moreover, 
techniques such as HREELS can be used to probe the hybridization of the substrate's phonon 
modes with graphene's ix plasmon under controlled doping [26, 28], and further explore the 
role of the nearby SPE continuum as a region in which such hybridized modes can undergo 
Landau damping. Metallic substrates, on the other hand, have attracted attention only recently 
as a possibly interesting dielectric environment for graphene [10, 11]. It has been recognized 
for many years that metal surfaces can exert strong dissipative and conservative (image) forces 
on moving external charges, especially in the grazing scattering geometry [31]. In particular, 
low-energy SPEs in a metal have been shown to give rise to strong friction on particles moving 
slowly over metal surfaces [20, 32, 33]. It is therefore tempting to explore how this concept 
of friction would be affected by the presence of graphene on a metal substrate, resulting in 
an overlap of the metal's SPEs and graphene's intra-band SPEs. Similar effects have been 
explored using a non-local description for a Cu overlayer on a Si substrate [34]. 

Consequently, we believe that studying both the loss function of the graphene-substrate 
system and the stopping force on an external moving charge can be used to explore the effects 
of substrate's dynamic response. There are many model parameters that may be relevant 
to the description of such processes. In our previous work [35], we have concentrated on 
amending the MDF-RPA dielectric function for graphene with a constant damping rate, 7, 
as well as on analyzing the effects of the equilibrium charge carrier density in graphene, 
n, and the graphene-substrate gap, h, for a substrate treated in the static mode. In this 
work, we consider further refinements to the MDF-RPA dielectric function for graphene, 
including a dispersing damping rate and the LFE correction. However, our main focus is 
on the inclusion of the dynamic and non-local properties of a substrate's dielectric response, 
for which we consider two specific examples of low-frequency excitations: a SiC substrate 
with a pronounced transverse optical (TO) phonon, and an Al substrate with a well-defined 
spectrum of low-energy SPEs. In doing so, we reduce the parameter space of our model by 
considering only two equilibrium charge densities: n = 0, characterizing intrinsic or undoped 
graphene, and n = n , where n = 10 13 cm~ 2 , as a value very close to graphene's charge 
carrier density for both the SiC [17] and Al [11] substrates. We also take h = 1 A as a 
reasonable estimate for the graphene-substrate gap [11, 35]. 

After outlining the basic theory in the following section, we present results comparing 
our model with the HREELS experiment of Liu et al. [17]. We then consider the effects of 
various model parameters on the stopping force for free graphene, followed by examples of 
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the stopping force for graphene on a SiC substrate and the friction coefficient for graphene on 
an Al substrate. A brief discussion of these results is provided in the concluding section. We 
use Gaussian electrostatic units throughout this paper unless otherwise stated. 



2. Theory 



We use a Cartesian coordinate system with coordinates {r, z}, where r = {x, y}, and assume 
that graphene is placed in the z = plane. A semi-infinite substrate with a bulk dielectric 
function e su b(k, uj) is assumed to occupy the region z < —h underneath graphene, while the 
region z > —h is assumed to be a vacuum or air [36]. The 3D wave vector describing the 
excitation modes in the substrate can be decomposed as k = {q, k z }, where q = {k x , k y } are 
the components parallel to the substrate surface and k z is the component perpendicular to the 
substrate surface. We further assume that an external charge with density p e xt( r , z, t) moves 
along a classical trajectory localized to the region z > above graphene. Following Ref.[36], 
we can express the induced potential $i n d(r, z, t) in the region above graphene by using the 
2D Fourier transform with respect to the surface components (r — > q) and with respect to 
time (t — > uj) as 



e(q,^J 



2tt 



S(q,u) 
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where e(q, uj) is the dielectric function of the combined graphene- substrate system and 



S(q,uj) 
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is the structure factor of the external charge distribution. The dielectric function of the system 
can be written as 

e(q,w) = e hg (q,u) + V c (q)U(q,uj), (3) 

where eb g (g, uj) is the effective background dielectric function of the substrate, Vc{q) = 
2%e 2 /q is the Fourier transform of the 2D Coulomb interaction, and H(q, uj) is the polarization 
function for free graphene. 

Assuming that the substrate is isotropic so that its dielectric function only depends on the 
magnitude of the wave vector, k = q 2 + fcf, one can use the SRM to handle the electrostatic 
boundary condition on the substrate's surface [29]. Combining the steps of the derivation 
outlined in Ref. [37] with those in [36], we obtain 



eb g (g,w) 



e s (q,uj) 
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where e s (q, uj) is the effective surface dielectric function of the substrate, defined by 

~q r° dk K 
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The expression (e s — l)/(e s + 1) in Eq. (4) represents a non-local surface response function 
for a substrate that exhibits strong dispersion, but in cases where the dispersion in the 
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substrate can be neglected so that e su b(A;, uj) ps e sub (^), one recovers the local approximation 

e s (q,u) « e suh (u). 

As an example of a non-local surface dielectric function, one may consider the low- 
frequency expansion of Lindhard's dielectric function, [20, 32, 33] 

tv^2 

e suh (k, uj) « 1 + + i^r H(2K F - k) , (6) 

where K^-p = \Z3Q p /Vf is the Thomas-Fermi (TF) inverse screening length of a 3D 
electron gas having a volume density N, f2 p = '4ire 2 N/m e is its bulk plasma frequency, 
Kp = m e Vp/h = (37i 2 N) 1 ^ 3 is its Fermi wavenumber, and m e is the mass of an electron. 
This dielectric function describes the low-energy SPEs of a metal subject to the cut-off 
k < 2Kp, as implied by the Heaviside function H (2Kp — k) in Eq. (6). We note that Eq. (6) 
has often been used in Eq. (5) to model various surface precesses on metals [20, 32, 33]. 

On the other hand, the local approximation may be suitable for cases such as a strongly 
polar insulating substrate characterized by non-dispersing TO phonon modes with frequencies 
co>TOi an d damping rates 7TOi> m which case one may use a dielectric function of the form 
[26, 27, 28] 

esubM = eoo + £ ~ r, (7) 

where = lim^oo e su b(c<j) is the high-frequency dielectric constant of the substrate and the 
oscillator strengths, fa, satisfy the relation ^\ fi = e — eoo, where e = e su b(0) is the static 
dielectric constant of the substrate [26, 27]. 

We note that in the literature on graphene it is common to assume a zero gap {h = 0) 
and to describe the substrate by only the static dielectric constant [2, 22, 23], for which Eq. 
(4) gives an effective background dielectric constant e^q, to) — > e e fr = (e + l)/2. In this 
case, a simple description of the screening of electron-electron interactions in graphene can 
be quantified by the Wigner-Seitz radius, r s = e 2 /(€ e shvp) [23], where vp ~ c/300 is the 
Fermi speed of graphene and c is the speed of light in vacuum. We note that the limiting case 
of free graphene (e^ g (q, uj) = 1) can be obtained by letting h — >• oo in Eq. (4) or by setting 
e = 1 in the zero gap case with a static substrate, but the focus of this work is on the effects 
of substrate including a finite gap, non-local effects, and its dynamic response. 

For the dielectric response of graphene in Eq. (3), we use the polarization function for 
graphene's n electron excitations in the MDF-RPA [22, 23], H(q, to), as well as the corrections 
to the polarization function for the LFE and the finite lifetime of the excitation modes of 
charge carriers in graphene. The LFE are described in the Hubbard approximation (HA) by 
replacing the MDF-RPA polarization function U(q, uj) with IIlfe^, = [1 — G(q)} n(g, to), 
where G(q) = (q/gd)/\/q 2 + kp and = 4 is the degeneracy factor of graphene's n 
electrons [2]. The finite lifetime of the excitation modes of charge carriers in graphene is 
treated by introducing a finite damping rate, 7, in the MDF-RPA polarization function through 
Mermin's procedure [24, 25], whereby one replaces H(q, to) with 

n M {q,w,V = 77 f tt 1 ... . . I — » (8) 
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where the static limit of the RPA polarization, H s (q) = lim^^o n(g, cu), is given elsewhere 
[22, 23]. 

In the case of a point charge Ze moving parallel to graphene with velocity v and at a 
fixed distance z > 0, we obtain S(q, a;) = 2irZe S(cu — q • v) e~ qz ° from Eq. (2). Using this 
expression in the inverse Fourier transform of Eq. (1), we can evaluate the stopping force on 
the point charge Ze from the definition 

F. = -Ze+-VQZ d {T,z,t)\^^ , (9) 
where v = v/v, giving [36] 



o 72„2 roc pqv 

F s = -^-\ dqe~ 2 ^ du> , U 9 



(10) 



71" v Jo Jo ^q 2 v 2 - u 2 |_ e (<?,^). 

We note that $X[—l/e(q,u)] represents the loss function of the system that, apart from a 
kinematic prefactor describing the scattering event, is directly accessible in the HREELS 
experiment [17, 38]. 

Finally, for cases in which the loss function in Eq. (10) can be expanded to the first 
order in frequency as $$[—l/e(q,u)] ~ ojT{q), one can define a friction coefficient r) for a 
sufficiently slow particle moving parallel to graphene through the equation F s « —rjv. In 
such cases, the friction coefficient is given by 

1 f°° 

7]=-Z 2 e 2 dqq 2 e- 2qzo F{q). (11) 



2 

For example, this procedure has been performed using Eqs. (5) and (6) to obtain the friction 
coefficient for slow particles on metal surfaces [20, 32]. A similar approach can also be used 
for graphene in the limit of vanishing damping rate, giving 

U(q, u) « II s (g) + H (2k F - q) J^-^-l, (12) 

where k F — y/Tm is graphene's Fermi wavenumber. A different but analogous expression can 
be obtained from Eq. (8) in the case of a finite damping rate [35]. 



3. Results 



3.1. Comparison with HREELS experiment 

It has been known for some time that the damping rate 7 of elementary excitations in a 
many-electron system can strongly influence the friction forces on slow particles [20, 32]. 
However, the exact value of the damping rate is generally difficult to measure, and in the case 
of graphene is presently unknown. Nevertheless, 7 is often used as an empirical parameter 
in the RPA response function modified by Mermin's procedure [24, 25]. Therefore, to obtain 
a reasonable estimate for the value of 7 for graphene, we compare the experimental data of 
Liu et al. [17] for the HREELS spectra of doped graphene on an SiC substrate with the loss 
function S$[—l/e(q, u)] obtained from Eq. (3) by using the Mermin polarization function for 
graphene, which is given in Eq. (8). In addition to constant values of 7, we also explore a 
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damping rate with a linear dispersion of the form 7 = v c q, where v c is a constant speed. 
The gap height h is also treated as an empirical parameter with a reasonable value of 1 A. 
We note that our comparison with the HREELS experiment is only tentative because we do 
not include the low-frequency contributions from scattering kinematics or temperature effects 
[38]. Nevertheless, we expect that the loss function obtained from the Mermin polarization 
function with a suitable choice of 7 will give the correct order of magnitude for spectral widths 
at frequencies that are not too low. 

The effects of a SiC substrate on graphene are not as well understood as those of a Si0 2 
substrate, but without entering the current debate [8, 9] we simply neglect any changes in 
graphene's 7r-band structure that may result from a hybridization of its ix orbitals with the 
substrate. In addition to treating the SiC substrate in the static mode with dielectric constant 
eo = 9.7, we reproduce some of the low-frequency features from the HREELS spectra [17] 
in a qualitative manner by including the dominant TO phonon mode of SiC in Eq. (7) with 
eigen-frequency co> TO = 97 meV [27], damping rate 7to =10 m eV [39], and high-frequency 
dielectric constant = 6.5 [27]. Also, by including the TO phonon in our comparison with 
the HREELS data, we gain some insight into the dynamic response of an insulating substrate 
prior to using it in calculations of the stopping force on a moving charge. 

In Fig. 1, we display a comparison of the HREELS data [17] for wavenumbers ranging 
from 0.008 to 0.126 A 1 with the loss function obtained from the MDF-RPA polarization 
function with 7 = (solid lines) and the loss functions obtained from the Mermin polarization 
function with 7 = v F k (dashed lines) and 7 = v F q (dotted lines), where k = v /vf% and 
n = 10 13 cm~ 2 . We also show cases for which the SiC phonon is included in the model (thick 
lines) and for which the SiC substrate is treated in the static mode (thin lines). The equilibrium 
charge carrier density in graphene is set at n = 1.9 x 10 13 cm -2 (hence k F ~ 0.08 A" 1 , giving 
the Fermi energy e F = hv F k F ~ 570 meV) to match experimental conditions [17]. Note that 
since the HREELS data is scaled to arbitrary units of intensity, the loss functions with finite 7 
are scaled so that their maximum peak heights coincide with those from the experiment. 

From Fig. 1, it can be seen that the best fit to the data can be obtained with the constant 
value 7 = v F k , while the best fit of the linear damping rate, 7 = v F q, tends to underestimate 
the widths of the experimental spectra at long wavelengths. Moreover, it appears that the 
inclusion of the substrate's TO phonon hints at the experimentally observed non-dispersing 
feature at approximately 100 meV. For the 7 = case, the two sets of narrow peaks indicate 
the position of graphene's it plasmon at long wavelengths before it crosses into the inter- 
band SPEs continuum at q « 0.07 A -1 and becomes broadened from Landau damping. 
Interestingly, the plasmon peak is broader in the case with the substrate phonon than in the 
case of a static substrate. This is a consequence of the hybridization of the substrate phonon 
with damping rate 7x0 = 10 meV and graphene's n plasmon with a vanishing damping rate. 
Therefore, using the known value of 7 TO = 10 meV as a "marker", one can surmise that the 
peaks shown in Fig. 1 with thick solid lines at long wavelengths are, in fact, the substrate 
phonon that is promoted due to an avoided crossing with the n plasmon [16]. This idea will 
be discussed in more detail shortly. 

Since the HREELS experiment [17] has been designed to determine the plasmon 
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dispersion in graphene on a SiC substrate, in Fig. 2 we show the peak positions deduced 
from the experimental spectra and from the various theoretical models displayed in Fig. 1. 
In addition, we show the peak positions resulting from the polarization function for graphene 
with the LFE correction, n L FE(<7? w), m the 7 = case (dash-dot lines) [2]. While it is clear 
that the LFE correction is unable to reproduce the widths of the experimental spectra, Fig. 2 
also demonstrates that its effects on the plasmon dispersion are hardly discernible from the 
case of the MDF-RPA polarization with 7 = 0. Fig. 2 further confirms that the case 7 = vpk 
provides the best fit to the experimental plasmon dispersion. For the static substrate case (thin 
lines), one notices that the peak positions approach the value c<j = 0asg^0ina manner that 
depends on damping rate, so that a typical 2D plasmon with uj ~ y/q is produced for 7 = 
and 7 = v c q, while a constant 7 appears to give rise to a quasi-acoustic relation uj ~ q when 
q -> 0. 

Finally, in all models of the plasmon peaks, the inclusion of the substrate phonon (thick 
lines) appears to capture the tendency of the experimental data to approach a finite frequency 
rather than vanish as q — > 0. A closer inspection shows that the frequency in question is near 
the value of 116 meV, which corresponds to the Fuchs-Kliever (FK) surface phonon at the 
frequency ujfk = ^to \/( e o + 1) / (^oo + 1) [27, 39]. This behavior is commensurate with 
the concept of a phonon-plasmon hybridization [26], which is illustrated in Fig. 2 by two 
thick solid lines showing the zeroes of the dielectric function for the MDF-RPA polarization 
function with 7 = 0. These two lines display an avoided crossing near the point where the 
thin solid line (representing graphene's n plasmon dispersion for a static substrate) would 
cross the FK phonon frequency ujfk ~ 116 meV. One can see from Fig. 2 that the substrate 
phonon frequency starts at uj fk for q = and is promoted to higher frequencies as q 
increases, following a dispersion curve similar to that of graphene's n plasmon in the case 
of a static substrate. This promoted phonon dispersion enters the inter-band SPE continuum 
at uj > v F (2k F — q) with u > vpq. On the other hand, graphene's plasmon starts from 
the origin of the (q, uj) plane and initially follows the standard long wavelength plasmon 
dispersion 0J p {q) = v F \J2r s kpq [22, 23, 36], but quickly levels to the TO phonon frequency 
co> TO = 97 meV as q increases. As a consequence of this avoided crossing [16], the plasmon 
branch is pushed into the intra-band SPEs continuum of graphene at uj < Vpq, where it 
broadens after crossing the point q c « ujto/ v f ~ 0.015 A -1 from Landau damping. We 
note that this complicated scenario is depicted in Fig. 2 for n = 1.9 x 10 13 cm" 2 , which 
is the spontaneous doping density of graphene on an SiC substrate [17], but by changing 
graphene's charge carrier density with a gate potential, for example, one can alter the picture 
of the plasmon-phonon hybridization quite significantly. 

3.2. Stopping force with a dynamic substrate 

After the discussion on the low-frequency excitations in graphene on an insulating substrate 
with a prominent TO phonon, we now explore how the various model parameters of this 
system affect the stopping force on a point charge, defined in Eq. (10). Given the complexity 
of the system, we first analyze the case of free graphene by letting h — > 00. In Fig. 3, we 
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show the velocity dependence of the stopping force on a proton (Z = 1) moving a distance 
z = 10 A away from graphene for the cases of the MDF-RPA polarization function with 
7 = (solid lines), the Mermin polarization function with 7 = uf^o (dashed lines) and 
7 = v F q (dotted lines), and the LFE polarization function with 7 = (dash-dot lines), where 
ko = v^™o an d ^0 = 10 13 cm" 2 . Results are provided for intrinsic graphene (n = 0) and for 
graphene with a charge carrier density n = n . From Fig. 3, one can see that all models give 
qualitatively similar results at high speeds, but at low speeds there are significant differences 
in the behaviour of the stopping force due to variations in the charge carrier density n and the 
damping rate 7. In particular, a threshold velocity at v = vp exists for the cases of intrinsic 
graphene with 7 = due to the lack of intra-band SPEs and the presence of inter-band SPEs 
in the lo > vpq region of the (q, uj) plane [22, 23]. Furthermore, all cases with 7 > or 
n > result in a linear dependence of the stopping force on the particle speed, which can be 
described using the concept of friction. The slopes of these curves do not appear to be affected 
significantly by the choice of a constant or linear 7, but the slopes do depend significantly 
on the density, as discussed previously [35]. It also appears that the LFE correction to the 
MDF-RPA polarization function has a small effect on the stopping force for speeds in the 
range v < vf- Therefore, we may conclude that the major influences on the polarization of 
graphene by an external charge are determined by the equilibrium charge carrier density n 
and the damping rate 7, which can be best described by a constant value. 

We now focus on analyzing the effects of the TO phonon in a SiC substrate on the 
stopping force of a proton moving a distance z = 10 A away from graphene with a graphene- 
substrate gap height of h = 1 A. Given our results for free graphene, in Fig. 4 we present 
the stopping force for the damping rates 7 = and 7 = vpk , as well as for the two charge 
carrier densities n = and n = n = 10 13 cm -2 . In all cases, the inclusion of the phonon 
(thick lines) results in a surprisingly large increase of the stopping force for speeds in the 
range v > 2v F when compared with the case of a static substrate (thin lines). This is due to 
the fact that as the speed increases, the integration range < uj < vq in Eq. (10) probes a 
larger portion of the promoted phonon dispersion curve shown in Fig. 2, which approaches 
a finite frequency rather than vanishing as q — > 0. However, we are more interested in the 
stopping force at low speeds (v < vp) for which the phonon-plasmon hybridization produces 
some peculiar effects. We first analyze the 7 = case in Fig. 4(a). For intrinsic graphene, 
which does not support a plasmon mode, one can see from the inset that the inclusion of the 
substrate phonon gives rise to a small but non-vanishing stopping force in the subthreshold 
region due entirely to the non-dispersing FK phonon at frequency ufk ~ 116 meV. For the 
finite density n = n , the inset in Fig. 4(a) shows that the linearly increasing stopping force in 
the static substrate case is augmented by a contribution from the original ix plasmon, which is 
pushed into the intra-band SPE region due to an avoided crossing with the substrate phonon 
and traverses this region at the constant frequency c^to ~ 97 meV without dispersion but 
with broadening from Landau damping. It is important to note that the inclusion of the TO 
phonon destroys the linear velocity dependence of the stopping force at low speeds, making 
the concept of friction inappropriate for the 7 = case. Moving to the case of 7 = vpko, 
which is shown in Fig. 4(b), one notices that the large peak in the magnitude of the stopping 
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force has been broadened and suppressed. Furthermore, although the effects of the substrate 
phonon are somewhat softened at low speeds, they still result in a non-linear dependence of 
the stopping force on the velocity for both n = and n = n . 

Finally, we analyze the stopping force of a slow charge moving parallel to a layer 
of graphene on an Al surface. From a first-principles calculation, it has been shown that 
graphene's ir electron band is not destroyed by a hybridization with aluminum, leaving its 
elementary excitation modes intact [11]. However, aluminum is close to an ideal metal with 
a volume density of quasi-free electrons given by N = 3/(47ri? 3 ), where R s « 2.07 Ob 
and a B is the Bohr radius, and supports a high-frequency bulk plasmon as well as a well- 
defined continuum of low-energy SPEs. Without addressing the problem of a hybridization 
of graphene's n plasmon and aluminum's surface plasmon, we explore the interplay of the 
intra-band SPEs in graphene and the contribution of the SPEs in aluminum to the non-local 
surface dielectric function, given in Eq. (5). Since a linear expansion of both the MDF-RPA 
polarization function for graphene, given in Eq. (12), and the surface dielectric function for 
aluminum, obtained from Eqs. (5) and (6), is possible for low frequencies, we can apply the 
concept of the friction coefficient 7] defined in Eq. (11). We note that the q integral in Eq. (1 1) 
takes into account two wavenumber cut-offs: q < 2k F = 2^flva for graphene, and q < 2K F « 
3.5 A -1 for aluminum. Since typical charge carrier densities in graphene have k F <C K F , it 
is interesting to see how significant graphene's screening ability is in a combined system with 
a metal. Following Ref.fll], we take equilibrium density of graphene on aluminum to be 
n = n = 10 13 cm -2 , and the gap height to be h = 1 A. 

In Fig. 5, we show the dependence of the friction coefficient on the distance z for a 
proton moving above free graphene, a free aluminum surface (which is placed in the plane 
zq = — 1 A), and a graphene on aluminum system. We also show the friction coefficient 
for a particle moving above free graphene with n — 0. Although we only consider the case 
of vanishing damping, we note that the effects of finite damping can be quite large in both 
systems [20, 32, 35]. For the finite density n = no, one can see from Fig. 5 that while the 
friction from the free Al surface dominates at short distances, the friction from free graphene 
has a longer range and dominates for z > 4 A. Furthermore, the friction coefficient of the 
combined system is suppressed when compared to both the case of free Al (at short distances) 
and the case of free graphene (at long distances). In the combined system with intrinsic 
graphene, for which there is no direct contribution from graphene's intra-band SPEs to the 
friction, one can see that graphene's static screening reduces the friction from low-energy 
SPEs in aluminum. Although we do not explore the effects of the density any further, from 
the n = and n = n cases one can surmise that the friction of a combined graphene on 
metal system can be controlled by changing the doping level of graphene externally. 

4. Concluding remarks 

We have provided a simple theoretical model that describes the dielectric coupling of 
graphene, represented by a 2D gas of massless Dirac's fermions, and a semi-infinite 3D 
substrate, represented by a surface response function in a non-local formulation. The 
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underlying assumption that the electronic band structures of these two systems is not altered 
by their proximity has allowed us to study the hybridization of their excitation modes due 
to Coulomb interaction. An emphasis was placed on the role of the dynamic response of a 
substrate in the low-frequency excitations of the combined graphene-substrate system, which 
give rise to the energy loss of slow moving charges above graphene. A comparison of the loss 
function for the graphene on an insulator system with a HREELS experiment was used as a 
motivation to explore methods of improving the polarization function for free graphene in the 
RPA. After identifying the damping rate and the equilibrium density of charge carriers as the 
two main parameters in graphene's polarization function, we calculated the stopping force on 
a point charge moving parallel to graphene supported by a SiC substrate, which was described 
by an empirical dielectric function that included an optical phonon excitation mode without 
dispersion. As a result, strong effects from the hybridization of graphene's ir plasmon and 
the substrate's surface phonon were found on the low-velocity stopping force. As another 
example, a friction coefficient was calculated for a slow charge moving above graphene 
on a metallic substrate by using Lindhard's dielectric function to obtain a surface response 
function representing the low-energy single-particle excitations in aluminum in a non-local 
formulation. In the limit of vanishing damping rates for both graphene and aluminum, an 
interesting interplay was found between the low-energy single-particle excitations of both 
systems in the friction coefficient. In all cases studied, the low-energy excitation modes in 
graphene were found to be strongly coupled with those in the substrate, giving rise to various 
forms of the stopping force on external charges that can be of interest to studies such as 
the chemical reactivity of graphene. Particularly intriguing is the possibility of affecting the 
dynamic response of the combined graphene-substrate system by changing the charge density 
in graphene by external means, thereby enabling a control of friction at the nano-scale. 
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Figure 1. The MDF-RPA (7 = 0) and Mermin loss functions (7 7^ 0) versus the energy loss 
for graphene with a charge-carrier density n = 1.9 X 10 13 cm -2 supported on a SiC substrate 
with a gap height h = lA. Model results for the Mermin loss function include the best-fit 
constant and linearly dispersing damping rates 7/(^f&o) = 1 an <3 7/ = 1, respectively, 
where ko = y/nno and n = 10 13 cm~ 2 . Thick and thin lines show model results with 
and without the inclusion of the substrate's TO phonon, respectively, while symbols show the 
experimental HREELS data from Ref. [17]. 
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Figure 2. The peak positions of the MDF-RPA loss function (7 = 0), the Mermin loss function 
(7 7^ 0), and the LFE-corrected MDF-RPA loss function (7 = 0, LFE) shown as a function 
of the wavenumber q/kp for graphene with a charge-carrier density n = 1.9 x 10 13 cm -2 
supported on a SiC substrate with a gap height h = 1 A. Model results for the Mermin loss 
function include the best-fit constant and linearly dispersing damping rates "f/(vpko) = 1 and 
7/(wjr?) = 1, respectively, where fco = y/irno and no = 10 13 cm~ 2 . Thick and thin lines 
show model results with and without the inclusion of the substrate's TO phonon, respectively, 
while filled circles show the experimental HREELS data from Ref. [17]. 
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Figure 3. The stopping force from the MDF-RPA loss function (7 = 0), the Mermin loss 
function (7 7^ 0), and the LFE-corrected MDF-RPA loss function (7 = 0, LFE) shown as a 
function of the reduced speed v/vpofa proton (Z = 1) moving at a distance zo = 10 A above 
free graphene (h = 00). Results for the Mermin loss function include the best-fit constant 
and linearly dispersing damping rates 7/(^^0) = 1 an d l/( v F<l) = 1, respectively, where 
fco = v/t^o an d "0 = 10 13 cm~ 2 . Thin and thick lines represent the case of intrinsic graphene 
and graphene with a reduced charge carrier density n/n n = 1, respectively. 
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Figure 4. The stopping force from (a) the MDF-RPA loss function (7 = 0) and (b) the 
Mermin loss function with j/{vpk ) = 1, where k = ^nno and n = 10 13 cm~ 2 , shown 
as a function of the reduced speed v/vf of a proton (Z = 1) moving at a distance zq = 10 A 
above graphene supported on a SiC substrate with a gap height h = 1 A. Results are shown for 
intrinsic graphene (solid lines) and graphene with a reduced charge carrier density n/n® = 1 
(dashed lines), as well as for the model with (thick lines) and without (thin lines) the inclusion 
of the substrate's TO phonon. 
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Figure 5. The friction coefficient 77 shown as a function of the distance zq of a proton (Z = 1) 
moving above graphene on an aluminum substrate (solid lines), free graphene (dashed lines), 
and a free aluminum substrate (dotted lines). Where appropriate, thin and thick lines represent 
the case of intrinsic graphene and graphene with a reduced charge carrier density n/riQ = 1, 
where no = 10 13 cm~ 2 , respectively. 



